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Abstract. In ORFEUS II spectra of the sdO star BD +39 3226 
interstellar hydrogen and deuterium is detected. From Ly a 
profile fitting and a curve of growth analysis of the Lyman se- 
ries of H I and D I we derive the column densities Ng. = 
1.20±gJSj-10 20 cmT 2 and iV D = 1.45±g jg-lO 15 cm" 2 . From 
the analysis of metal absorption lines in ORFEUS and IUE 
spectra we obtain column densities for 1 1 elements. In addi- 
tion, we examine absorption lines of H2 for rotational exci- 
tation states up to J = 7. We find an H2 ortho-to-para ra- 
tio of 2.5, the fractional abundance of molecular hydrogen 
has a low value of log/ = —4.08 for a total amount of 
jV(H 2 ) = 4.8± 2 ;g-10 15 cm- 2 . The column densities of the ex- 
citation states reveal a moderate Boltzmann excitation temper- 
ature of 130 K and an equivalent excitation temperature for the 
excited upper states due to UV pumping of < 1800 K. 

Key words: Space vehicles - Stars: individual: BD +39 3226 
- ISM: abundances - ISM: molecules - Cosmology: miscella- 
neous - Ultraviolet: ISM 



1. Introduction 

The deuterium abundance in the Galaxy has been of interest 
since the first models of primordial nucleosynthesis were de- 
veloped. While it is unlikely that major amounts of deuterium 
could be produced after the Big Bang, it is sure that deuterium 
is almost completely destroyed in stars. This should lead to a 
decrease of its abundance in the progress of galactic evolution. 
Every measurement of D/H thus should give a lower limit 
for the primordial abundance and, according to nucleosynthe- 
sis models, an upper limit for the baryonic density in the uni- 
verse. Today, these limits are set by me asurem ents in quasar 
spectra. For example Levshakov et al. ( 1999 ) suggest a uni- 
form D/H ~ 4-10~ 5 for three lines of sight at redshifts of z 
between 0.7 and 3.6, but Molaro et al. ( |1999| ) find 1.5 • 10" 5 
from a line of sight with z — 3.5. Nevertheless knowledge of 
abundances in the galactic ISM can be important for tracing the 
evolution of the Milky Way. 



Since the Copernicus satellite in the 1970s no instrument 
was capable of high resolution spectroscopy in the 900 to 
1200 A range. There have been D/H measurements with the 
IUE and the HST-GHRS, but only the Ly a line was observable 
leading to a restriction to lines of sight with low hydro gen co l- 
umn densities. For an overview see e.g. Lemoine et al. ( 1999 ). 

The ASTRO-SPAS space shuttle platform housed 3 spec- 
trographs operating between 900 and 1200 A. Of those, IMAPS 
has an echelle spectrograph working between 930 and 1150 A 
with a resolution of A /AA ~ 1.5 • 10 5 . Jenkins et al. ( |1999| ) 
and Sonneborn et al. ( 1999| l performed D/H measurements 
with that instrument. Two spectrographs were attached to the 
ORFEUS telescope. The Berkeley spectrograph is designed for 
intermediate-resolution spectra in the range of 390 to 1220 A. 
The Heidelberg-Tubingen echelle spectrograph, equipped with 
a microchannel plate detector, gives spectra from 912 to 1410 A 
and allows investigations of the entire hydrogen and deuterium 
Lyman series with a resolution of A/AA ~ 10 4 (Kramer et al. 



1990). We use spectra of the latter to investigate deuterium. 
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The large number of rotational transitions of molecular hy- 
drogen found in the wavelength range from 1200 A up to the 
Lyman edge can be used for studies on the molecular gas along 
the line of sight. The distribution of molecular hydrogen in the 
local interstellar medium especially at higher galactic latitudes 
is still only rudimentarily known, since the majority of the 
Copernicus measurements pointed towards luminous targets in 
the galactic plane. From the H2 transitions column densities of 
the different H2 rotational excitation states up to J = 7 can be 
determined. These can be used to derive physical parameters 
as the excitation temperature and the ortho-to-para ratio of the 
molecular hydrogen. 



2. Observations and data reduction 

BD +39 3226 is a sdO star displ aying a pure helium line spec- 
trum in the optical (Giddings 198C ) in a distance of about 
270 pc[] at the galactic coordinates I = 65.00, b = +28.77. 
Its high radial velocity of ~ —279 km s _1 makes this star suit- 
able for studies of the local ISM, since stellar and interstellar 
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absorption lines are well separated. We analyse echelle spectra 
obtained during the ORFEUS //mission in Nov./Dec. 1996 and 
NEWSIPS reduced long- and short-wavelength high dispersion 
spectra taken from the IUE Final Archive (IUEFA)K The total 
wavelength range thus covered is about 900 to 3200 A. 

The ORFEUS spectrum was obtained in 4 pointings with 
a total of 5000 s integration time. The basic reduction of the 
data was performed by the ORFEUS team in Tubingen. Details 
about the ORFEUS instrumentation and the basic data handling 



and calibration are given in Barnstedt et al. (1999). We repeat 
the main features here. 

The wavelength calibration of the ORFEUS spectra is 
based on the interstellar spectrum in the ORFEUS target HD 
93521. The accuracy of this calibration is estimated as bet- 
ter than 0.05 A, but small systematic effects within the spec- 
trum cannot be excluded. The zero point of the wavelength 
calibration is based on the geocoronal Ly a emission. Since 
the ORFEUS aperture measures 20" in diameter, errors in the 
pointing may shift the wavelength zero point. Optical spectra 
of BD +39 3226 give the stellar heliocentric radial velo city as 
— 278. 7± 5 km s_Vbased on values by Dworetsky et al. ( 1982 ) 
and Giddings ( |1980| ) and one derived from observations with 
FOCES at the Calar Alto 2.2 m telescope. In the ORFEUS 
spectrum the mean radial velocity of 22 sharp stellar metal lines 
is -289.1 ± 1.5 kms" 1 . We corrected the ORFEUS spectrum 
for this shift. The new velocity zero point should be accurate 
within w 5 kms" 1 . 

The spectral resolution is intrinsically better than 10 4 , but 
pointing jitter (~ 2") and the addition of spectra from sev- 
eral pointings may cause some deterioration. We have filtered 
the spectra as provided by Tubingen with a wavelet algorithm 
(Fligge & Solanki 1997) effectively leading to a 3 pixel boxcar 
filter. We derive from our spectrum a resolution of A/ A A ~ 
0.9-10 4 . 

ORFEUS spectra are influenced by scattered light, which 
is implicitly corrected for by subtracting the intensities in the 
interorder space from the echelle order intensities. Near strong 
features residual effects may still be present (Barnstedt et al. 



1999). However, since the ORFEUS spectra are slightly tilted 



with respect to the detector grid, affected areas lie normally at 
some distance from such features. 

The resolution of the IUE spectra is also about 10 4 , which 
is equivalent to ~30 krns" 1 . While only one long wavelength 
spectrum (LWR 11789, exp. time 5040 s) can be found in the 
IUEFA, several short wavelength spectra are available (SWP 
15275, 48312, 48313, 48314, each with 3600 s exp. time). To 
improve the S/N ratio we added the four SWP spectra. The IUE 
spectra were smoothed with a 3 pixel boxcar filter. 

3. Column densities 

3.1. Data analysis 

We identified numerous interstellar absorption lines of H I, D I, 
different heavy elements and H2 . In the photospheric spectrum 
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Fig. 1. Profiles of metal absorption lines in the IUE LWR, 
IUE SWP, and the ORFEUS spectrum. The velocity compo- 
nents (heliocentric) at i> ra d — —75 kins" 1 (B) and i> ra d — 
—24 kms" 1 (A) are labeled 



only a few weak, sharp metal lines of C III, C IV, N III, N V, 
Si III, Si IV, P V, and S V can be identified besides the strong 
He II lines. The investigations on interstellar metal and deu- 
terium abundances make use of the standard curve of growth 
technique. Except for the case of the Ly a and i5 absorptions, 
the equivalent widths of the lines were measured either by a 
trapezium or by a Gaussian fit. The differences between the 
two methods are well below the typical errors in the equiva- 
lent widths. Uncertainties in the measurements occur because 
of noise and the choice of the continuum. The error due to the 
latter was estimated by determining the equivalent width for 
a lower and an upper limit of the continuum level, the photon 
statistics were taken into consideration by a formula adopted 
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from Jenkins et al. (1973 ). Both errors added quadratically give 
the uncertainties in W\ as quoted in Tables [j], ^, and |j. 

Some lines of S II, Si II, and N I, which lie between 1 190 
and 1390 A, were measured in both the IUE and the ORFEUS 
spectrum. The line profiles and absorption strengths turn out to 
be consistent. 

3.2. Metals 

In order to judge the velocity structure of the absorption it is 
reasonable to analyze the metal lines at first. We find 2 dif- 
ferent absorption components: one at a radial velocity (he- 
liocentric) of ~ — 24 kms" 1 (component A) and one at ~ 
— 75 kms" 1 (component B). The latter is rather weak, Sill 
equivalent widths are smaller than 24 mA with large uncertain- 
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ties, so an examination of component B can give only very un- 
certain results. Even if this component had a very low 6-value, 
the upper limits for the S II and Si II column densities would be 
roughly 1.6T0 14 cm~ 2 and 7.9-10 14 cm~ 2 respectively corre- 
sponding to a hydrogen column density of < 10 19 cm -2 . 

In the following we will concentrate on the cloud at 
—24 kms -1 . The small LSR velocity of this component sug- 
gests a local origin while component B at —75 kms -1 repre- 
sents most likely a cloud at larger distance. 

Table [j] lists the measured equivalent widths for compo- 
nent A. The lines of Fe II, Si II, and N I were used to determine 
the b- value of 5±1 kms -1 of the curve of growth for metals. 
Then column densities of other ions were obtained by fitting 
their equivalent width to the curve. Fig. || shows the curve of 
growth, Table || gives the resulting column densities. The er- 
rors take the uncertainties in the fe-value and in the individual 
equivalent widths into account. For ions with data points only 
in the flat part of the curve of growth the column densities have 
large uncertainties. 



3.3. Neutral hydrogen 

We determined the H I column density in two different ways. 

First we compared theoretical Voigt profiles convolved with 
a gaussian instrumental profile to the Ly a line in the ORFEUS 
and the IUE spectrum. This line is always fully damped, there- 
fore the 6-value is unimportant and in case of a single veloc- 
ity component only the column density remains as a param- 
eter. Even small changes in the column density have a clear 
effect on the profile. We estimate the accuracy in \ogN as 
±0.1. Problems arise because of the weaker component B at 
AA~ - 0.21 A and stellar Hell absorption at AA~ - 1.57 A 
which both are unresolved. Component B should have only a 
weak influence on the column density, probably < 0.04 dex. 
The stellar line was calculated from an atmospheric model 
(T cff = 45000 K, log 5 = 5.5, n(He)= 99%, n(H)= 1%) 
and included in the fit. The uncertainty in A^(Hl) due to the 
stellar model is small, because significant errors in the strength 
of the calculated stellar line would have made the fit profile 
asymmetric with regard to the measured profile. An additional 
background substraction was applied in both spectra of the Ly- 
man a line to correct for some residual intensity (w 5.5% in 
the ORFEUS spectrum) near the centre. The result is plotted 
in Fig. [|. It is only possible to give a total column density (for 
component A and B), which is N H = l.^l^-lO 20 cm" 2 . 

To confirm this value, we also applied the curve-of-growth 
analysis to the Lyman series from Ly 5 to Ly k, except for Ly 9 
and Ly l which seem distorted. The Ly (3 and Ly 7 lines have 
strong damping wings which, together with the further struc- 
ture of the spectrum, do not allow the determination of reli- 
able equivalent widths. Higher Lyman lines are also visible but 
may be blended with stellar He II, because the distance between 
these lines is smaller than the stellar radial velocity. Besides, 
near the Lyman edge it is not possible to set the continuum 
properly. We measured the equivalent width of Ly i5 by compar- 
ing the line with computed two-component-profiles (Voigt pro- 
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Fig. 4. Theoretical fits to the Lyman a line in the IUE and the 
ORFEUS spectrum. In both spectra the residual intensity near 
the centre of the absorption has been set to zero. The contribu- 
tion of the stellar Hell line located ~ 1.57 A towards shorter 
wavelengths as calculated from an atmospheric model is in- 
cluded in the fits and plotted seperately. The two shown fits 
represent log A r (Hl)= 20.0 and 20.1. Near the centre the geo- 
coronal emission peak is visible 



files convolved with the instrumental Gauss profile), one com- 
ponent for the H I and one for the D I line. For the other lines 
the damping wings are neglegible and we used two-component 
Gauss fits. 

We note that for the higher Lyman series lines (6-k) the in- 
strumental profile degrades the true absorption such that resid- 
ual light is expected near the bottom of the profiles (see Fig. ^ 
and 0). A calculation shows this to be at the level of up to 10%. 
In addition, also in these lines geocoronal emission is present 
but not readily recognizable in the profiles considered (in Ly- 
man a and (3 it is clearly present). Since the absolute level of 
the contamination is not reliably known we have refrained from 
corrections. 

In case of Ly e the Fe II line situated between the H I and 
the D I line at 937.652 A was modeled and used as a third, 
fixed component in the fit. An analogous procedure was neces- 
sary for the H 2 Werner Q(l), 4-0 line lying at A = 930.574A 



4 



H. Bluhm et al.: Deuterium and molecular hydrogen in the ISM towards BD +39 3226 




Fig. 2. The curve of growth for ab- 
sorption by metals on the line of sight 
to BD +39 3226 (component A) is 
shown. Plotted are only Si II, Fe II, and 
Nl which actually define the curve of 
growth. Filled symbols represent OR- 
FEUS data, open symbols those from 
the IUE. The drawn curves represent 
single cloud absorption with indicated 
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Fig. 3. The curve of growth for H I and 
D I towards BD +39 3226. In the damp- 
ing part the curves split, the lower curve 
is for Ly k, the upper one is for Ly S. 
Note that the D I points lie essentially 
on the Doppler part of the curve of 
growth 
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Fig. 5. Curve of growth for H I in detail. The data points and 
the corresponding theoretical curves are labeled 



between the H I and D I Lyman £ lines. Attempts to include 
velocity component B by additional fit components showed 
that it has neglegible influence on the line shape. Examples 
for the fits are shown in Fig.^. Though the velocity structure 
is not resolved in the Lyman series, a single-cloud curve of 



growth should be sufficient since there is one clearly dominant 
cloud. As expected the H I curve of growth has a significantly 
higher Doppler-velocity than the metals' curve due to the much 
smaller atomic mass of hydrogen. For a given b each measured 
equivalent width and its error correspond to a column density 
with an error depending on the slope of the curve of growth. 
We calculated log Ah for different 6-values from the weighted 
mean of the 5 column densities resulting from the 5 measured 
equivalent widths. The least mean square deviation is found 
for 6 = 10.5 ± 1 km s -1 , leading to Ah = 1.6±g ;|-10 20 cm" 2 . 
Fig. |^ and |] show the curve for H I. 

The Lyman a fit and the curve of growth analysis give con- 
sistent results, so a mean value of Ah = I^O+q^-IO 20 cm~ 2 
can be derived. 

3.4. Deuterium 

The absorption of deuterium was investigated in 5 lines (see 
Table ^]). The D I Lyman 8 and rj lines are shown in Fig. [7| 
Along with the H I data, the D I equivalent widths are plotted 
in Fig. [3| For D I Ly n only an upper limit can be given which 
is rather high due to the continuum uncertainty. The Doppler 
parameter is only of minor importance for the determination of 
the D I column density because the data points lie mainly on the 
linear part of the curve of growth. The deuterium data points 
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Fig. 6. The Lyman lines used in our curve of growth analysis in 
the ORFEUS spectrum. The positions of D and H absorption 
by component A (u ra d — —24) are marked 



seem to suggest a somewhat higher b- value than 10.5 km s -1 
but we used the same curve as for hydrogen since we expect 
bu < &h- The weighted mean of the column densities derived 
for b — 10.5 kms -1 from the four measured W\ values is 
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3.5. Molecular hydrogen 

The ORFEUS spectra also contain a large number of absorption 
lines from molecular hydrogen. Only one component is visible 
here. The average radial velocity of 18 lines measured in the 
echelle orders 50-59 is "Uheiio = —27.1 ± 0.9 krns" 1 , slightly 
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Fig. 7. For two of the Lyman lines a 500 kms" 1 wide section 
of the spectrum is displayed. Top: Deuterium and hydrogen 
Lyman 5 line in the ORFEUS spectrum (data points) with a 
two-component fit (Voigt profile convolved with instrumental 
Gauss). Bottom: Deuterium and hydrogen Lyman r\ line (data 
points) with a two-component Gauss profile. In both plots the 
positions of the velocity components A and B found in the 
metal lines are marked for the H I lines. Component B has only 
neglegible effect on the absorption profile 

different from the velocity found for the metal absorption lines. 
The metal radial velocity of Uhciio = —23.9 ± 1.0 kms -1 , was 
measured as the average value of 20 lines in the echelle orders 
42-60. Since we have not found any obvious systematic veloc- 
ity shift between different orders, a possible explanation may 
be the presence of an additional weak unresolved absorption 
component in the metal lines. 

We have determined equivalent widths using trapezium fits. 
No better quality of the results would be achievable from fit- 
ting gaussian profiles since most of the lines are only weak and 
do not show clear profiles. The equivalent width for rotational 
states J > 5 are similar to the strength of noise peaks, so only 
upper limits are determined. Results are presented in Table Q 

Column densities for the different rotational excitation lev- 
els are derived from these equivalent widths by fitting to a the- 
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oretical curve of growth (Fig. ||). The fe-value is restricted by 
the lower J levels to 2.5 to 3 kms -1 . The column densities 
for J = and 1 lie in the flat part of the curve of growth and 
therefore are sensitive to variations of the 6-value, enlarging the 
errors for their column densities. The higher rotational levels 
are located on the doppler part of the curve which in principle 
allows a good definition of column densities. However, most 
of these lines are weak and have larger errors in the equivalent 
widths which again leads to larger errors in the column densi- 
ties. For J > 5 only upper limits can be given. The resulting 
column densities are listed in Table 

The population of the lower rotational states of molecular 
hydrogen is determined by collisional excitation, following a 
Boltzmann distribution. The excitation temperature can there- 
fore be derived as 



T — 



E n — E„ 



2.303 k log (N m /g m ) - log (N n /g n ) ' 



(1) 



where the statistical weight gj is equal to (2 J + 1), multiplied 
by a factor of 3 for J odd in regard to the triplet nature of 
ortho-H2. For the upper levels UV photons have to be consid- 
ered as the primary source of excitation (Spitzer & Zweibel 



1974; Spitzer et al. 1974). An equivalent excitation tempera- 



ture can be derived here, but this has by no means the physical 
implications of an actual kinetic temperature. 

We have determined excitation temperatures by using an er- 
ror weighted least square fit. A temperature of T02 — 130 ± 1 K 
is derived for the lower levels J = to 2, for the upper levels 
we find an upper limit in the equivalent excitation temperature 
of T 37 < 1800 K. 

Looking at Fig. |l0[ where the column density Nj, weighted 
by the statistical weight gj, is plotted against the excitation en- 
ergy of the rotational state J, a clear distinction is visible be- 
tween the collisionally Boltzmann excited levels J < 2 and 
the UV photon excited levels J > 3. Although the upper levels 
scatter around the fit, a much smaller slope is obvious, com- 
pared with the lower levels. 

The H2 ortho-to-para ratio (OPR) derived from the by far 
most prominent rotational states J = and 1 is 7V(l)/iV(0) ~ 
2.5. The excitation temperature calculated using the lowest 
para-states J = and 2 is T02 = 130 K. This value essen- 
tially coincides with the temperature determined from the OPR 
of Tqi = 134 K. We therefore can assume that we see thermal- 
ized, purely collisionally excited gas in the excitation levels up 
to J = 2, f ollowing the Boltzmann distribution (see Dalgarno 
et al. "~ 



1973) 



For the higher excitation levels information is less clear due 
to the scatter in the data points. Fitting the ortho- and para- 
states seperately we get T 35 < 870 K and T i6 < 1400 K, so 
the excitation temperature for the upper states probably can be 
restricted to < 1400 K rather than < 1800 K as stated above. 
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Fig. 8. Overview of the spectra for some H2 transitions. The 
spectra shown are filtered by a wavelet algorithm. The H2 ab- 
sorption at Whoiio — —27 kms -1 is clearly visible for the tran- 
sitions with J < 4. At higher rotational states the equivalent 
widths are similar to the strength of noise peaks. Additional 
absorption lines marked are: a) Ly R(l), 2-0; b) Ly R(0), 4-0; 
c) Ly P(l), 4-0; s) features of probably stellar origin, which 
could not be unambigiously identified 



4. Abundances 
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4.1. Metals, dust, and H2 

The metal abundances in component A are given in Table ||. 
Since only a total hydrogen column density for components A 
and B could be measured, the abundances may be too low by 
< 0.04 dex. In presence of the relatively large uncertainties this 
effect seems neglegible. 

As can be seen in Table || the depletion factors log d = 
log(N/Nn) — log(N/N-n)Q of the elements vary between 
—0.3 and —1.7. This moderate depletion indicates a low abun- 
dance of dust along the line of sight. 

The total amount of H2 seen on this line of sight is 



JV(H 2 ) 



-1.6 



10 cm , which gives an abundance rel- 



ative to H I of -4.40 



+0.17 

-0.20 



dex. With this low value the total 



amount of hydrogen, 7V(Hl + H 2 ) = [AT (Hi) + 27Y(H 2 )] = 
1.2-10 20 cm -2 , is hardly larger than the value given above for 
just Hi. 

The e xtincti on towards BD +39 3226 is given by Dworet- 
sky et al. ( 1982 ) as E(B — V") = 0.05, which again indicates a 
low dust abundance. In addition, the fractional abundance of 
molecular hydrogen / = 2A(H 2 )/ [TV (Hi) + 2/Y(H 2 )] has 
a low value of log / = —4.08. The gas to color excess ra- 
tio for our line of sight of N (H I + H 2 ) /E(B - V) = 2.3 • 
10 21 cm _2 mag _1 is slightly lower than the typical value for 
the galactic intercloud gas of 5.0 ■ 10 21 cm -2 mag -1 given by 



Bohlin et al. (1978). This deviation, though, can be explained 
by uncertainties in the E(B — V) value and should not be over- 
interpretated. We can conclude that most of the absorption seen 
on this line of sight arises from intercloud gas with only little 
contribution from the clouds having T ~ 130 K. 



Savage et al. (1977) have generated a plot of log/ vs. 
E(B — V) for their survey of 70 stars. For lines of sight with 
E(B — V) < 0.08, they find only very small fractions of H 2 
down to values log / of — 6. This is in line with the tight con- 
nection between the existence of dust grains and of molecules 
forming on their surfaces. Our low value fits well into this 
group of stars of low fractional abundance. 



In their plot, Savage et al. (1977) show theoretical curves 
from equilibrium models of Black for different hydrogen den- 
sities, radiation densities, and molecule formation rates. Ac- 
cording to this, our value is described by Black's "Model 2", 
which assumes a density n(Hl + H 2 ) of 100 cm~ 3 at a tem- 
perature of about 100 K, in accordance to the temperature de- 
rived above. The radiation density at 1000 A has then a value 
of 4 ■ 10 -17 erg cm -3 A -1 if the molecule formation rate is 
10" 17 crr^s" 1 . 



4.2. Deuterium 

From the H I and D I column densities we obtain a deuterium 
abundance on this line of sight of D/H = l-2±g:|-10 -5 . We 
also tried to identify absorption lines of HD, but without suc- 
cess. Taking the H 2 column density of only 4.8-10 15 cm~ 2 one 
would expect HD column densities as low as about 10 11 cm~ 2 , 
which is definitely below the detection limit of ~ 10 13 cm -2 
for the spectra available. It can be ruled out that any significant 
amount of deuterium is hidden in HD. Deuterium depletion by 
interactions with dust grains is unlikely in a warm ISM compo- 
nent with low gas-to-dust ratio. 

The abundance is a value for the whole line of sight, possi- 
bly an average over more than one cloud with different D/H 
ratios. The uncertainties are rather large, but the value is reli- 
able within its limitations. There was no need for assumptions 
about the velocity structure of the absorption (as it is necessary 
for line modelling) since the hydrogen column density deter- 
mination is mainly based on fully damped absorption while the 
deuterium lines lie in the linear part of the curve of growth. In 
both cases the curve of growth is insensitive to effects arising 
from clouds of different 6-values and column densities. Fur- 
thermore, a model of the stellar spectrum was used only for 
Lyman a fitting. The uncertainty arising from modelling is ne- 
glegible because the centres of the stellar and interstellar ab- 
sorption lines are well separated due to the high radial velocity 
of the star. 
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Fig. 10. Excitation plot for the H2 column densities. Plotted 
is the statistically weighted column density vs. excitation en- 
ergy of the rotational states. A clear separation of collision- 
ally excited levels (J < 2) and of UV photon excited levels 
(J > 3) is visible. The lines denote least square fits to an as- 
sumed Boltzmann distribution for levels m to n with the ap- 
propriate temperature T mn . The error given for the Boltzmann 
excitation temperature X02 is the standard deviation of the fit. 
For T37 only an upper limit can be given 

Varying scattered light may affect the measured equivalent 
widths, as mentioned in Sect. 2. However, we do not find ev- 
idence for systematic background substraction errors beyond 
the 5% level in all analyzed H and D lines and the scatter of the 
equivalent widths is not larger than expected for the estimated 
errors. The effect of systematic errors would be rather small. 
For example, if all H and D equivalent widths were larger by 
10 % than measured due to wrong background substraction, 
the D I column density also would rise by « 10% while the 
H I column density resulting from curve of growth analysis and 
Lyman a fitting would change even less. Thus the D/H ratio 
would be increased by less than 10 % and would remain well 
within the given range. 

5. Concluding remarks 



McCullough (1992) finds all Copernicus and IUE ISM deu- 



terium abundance measurements to be consistent with D/H = 
(1.5 ± 0.2)-10~ 5 . This value was confirmed for the local in- 



terstellar cloud by GHRS observations towards a Aur (Lin- 
sky et al. |1995[ ) and HR 1099 (Piskunov et al. |1997| ) yielding 



D/H = UHT^-IO- 5 and 1.46 ± 0.09-lQ-^respectively 



For a summary of GHRS measurements see Linsky (1998) 



McCullough's list of included data contains 6 measurements 
towards nearby stars (d< 14 pc) observed only in Lya and 8 
Copernicus observations of each 2 to 5 Lyman lines towards 
more distant stars (d> 79 pc). The mean of the results of the 
latter observations is D/H ~ 1.4-10~ 5 . Our result based on 
5 deuterium lines lies below this value though it is consistent 
within its uncertainty. 

Recent dete rmina tions of D/H for example by Vidal- 
Madjar et al. (|l99S|) towards G191-B2 B usin g the GHRS 
(D/H = 1.12+^|-10- 5 ), Golz et al. ( |1998| ) in a prelim- 
inary analysis of ORFEUS II observatio ns of BD +284211 
(D/H = O.St^-lO- 5 ) or Jenkins et al. ( |1999| ) using IMAPS 



spectra of 8 Ori (D/H = 0.741^- 10~ 5 ) indicate smaller val- 
ues towards early type stars with d > 69 pc. As Vidal-Madjaret 
al. and Jenkins et al. state, a possible explanation may be blend- 
ing of the Local Interstellar Cloud (LIC) contribution with that 
of more distant deuterium-poor clouds. Small scale spatial in- 
homogeneities of D/H could be the result of the production 
of significa nt am ounts of deuterium in stellar flares, as Mullan 
& Linsky ( 1999 ) suggest. This discussion remains open since 
Sahu et al. ( |1999| ) found D/H = 1.60±g;|f • 10~ 5 towards 
G191-B2B from a STIS spectrum, a value consistent with a 
homogeneous D /H in the LISM. 

However this may be, the various determinations of D/H 
all are consistent with D/H w (1.2 ± 0.3)-10~ 5 . The uncer- 
tainty in each individual value is influenced by the spectral res- 
olution of the instrument, by the number of D lines used, and 
by the actual D column density. Our result from the new in- 
strument ORFEUS using 5 lines at 30kms _1 resolution adds 
to these determinations . 
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Table 1. Metal absorption line equivalent widths measured in 
component A (at w rad ~ -24 kms -1 ) with ORFEUS and IUE 



Ion 



Wavelength 

[A] 



log(/A) 



W A (ORFEUS) 
[mA] 



W A (IUE) 
[mA] 



Fell 


1096.8769 


1.545 


25.6 ± 8.0 




- 








1125.4477 


1.092 


12.8 ± 6.4 




- 








1144.9379 


2.080 


64.1 ± 6.9 




- 








1608.4511 


1.998* 


- 


112 


1 ± 


18 


8 




2344.214 


2.410 


- 


161 


9± 


16 


4 




2374.4612 


1.889* 


- 


101 


4 ± 


15 


8 




2382.039 


2.967 


- 


163 


9 ± 


15 


9 




2586.6500 


2.248* 


- 


161 


3 ± 


26 


7 




2600.1729 


2.765 


- 


189 


1 ± 


20 


8 


Sill 


989.8731 


2.119 


66.0 ± 15.5 




- 








1190.4158 


2.474 


80.2 ± 11.1 


79 


4± 


7 


9 




1193.2897 


2.775 


101.1 ± 10.0 


94 


5± 


7 


1 




1260.4221 a 


3.104 


98.8 ± 16.4 


92 


9± 


11 


4 




1304.3702 


2.050* 


69.6 ± 15.2 


70 


9± 


6 


5 




1526.7066 


2.225* 


- 


105 


4± 


8 


1 




1808.0126 


0.596* 


- 


33 


3± 


8 


4 


Sim 


1206.500 


3.311 


69.6 ± 10.9 


64 


0± 


12 


9 


Nl 


964.6256 


1.059 


31.2 ± 8.7 




— 








1134.9803 


1.660 


68.5 ± 10.2 




- 








1199.5496 


2.202 


115.2 ± 10.4 


80 


0± 


5 


2 




1200.2233 


2.026 


82.1 ± 7.6 


62 


5± 


7 


8 




1200.7098 


1.725 


54.9 ± 9.2 


70 


5± 


9 


6 


Ol 


924.952 


0.168 


49.4 ± 16.7 




- 








936.6295 


0.544 


67.2 ± 11.9 




- 








950.8846 


0.174 


57.7 ± 7.0 




- 








976.4481 


0.508 


65.1 ± 11.1 




- 








1039.2304 


0.980 


70.4 ± 9.7 




- 








1302.1685 


1.804 


b 




6 






CI 


1277.2454 


2.091 


- 


20 


5± 


3 


4 




1560.3092 


2.099 


- 


23 


4± 


8 


8 




1656.9283 


2.367 


- 


45 


9± 


7 


7 


Cn 


1036.3367 


2.124 


94.4 ± 10.7 












1334.5323 


2.235 


134.0 ± 15.6 


142 


4 ± 


4 





SlI 


1250.584 


0.834 


34.2 ± 10.0 


36 


3± 


5 


4 




1253.811 


1.135 


22.2 ± 7.5 


46 


4± 


6 


3 




1259.519 


1.311 


50.0 ± 10.9 


59 


2 ± 


8 





Mgl 


2852.9642 


3.718 




79 


9± 


17 


1 


Mgll 


2796.352 


3.235 




239 


7± 


14 


9 




2803.531 


2.934 




250 


6± 


20 


6 


Mnll 


2576.877 


2.956 




56 


7± 


11 


2 




2594.499 


2.847 




51 


9± 


13 





Nill 


1370.132 


2.254 




6 


9± 


5 







1741.549 


2.256 




14 


5± 


5 


8 


Znll 


2026.136 


2.996* 




43 


9± 


14 


2 




2062.664 


2.723* 




10 


3± 


6 


9 


Alii 


1670.7874 


3.486 




92 


0± 


6 


7 



Wavelengths and /-values are f rom Morton (1991), except for 
from Savage & Sembach ( 199C ) 



taken 



a Line blended with Fell at 1260.533 A. The equivalent width given 
here has already been corrected for a contribution of 31.5 ± 4.6 mA 
derived from the determined Fe II column density 
b The 1302 O I line is in ORFEUS spectra blended with geocoronal O I 
emission, in IUE spectra with a reseaux 
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Table 2. Metal column densities N [cm 2 ] and abundances. Table 4. H2 equivalent widths and column densities Nj [cm 2 ] 



Solar values are taken from de Boer et al. (1987) 



Element 


1 AT 

log IS 


log (A/ Ah) 


log(AyA'H)© 


C 


16.40 ±0.75 


-3.70 ±0.75 


-3.4 


N 
O 

Mg 

Al 

Si 


14.75 ±0.25 
16 40 +°- 75 
14.95 ± 0.80 
13.00 

14.80 ± 0.20 


-5.35 ±0.25 

°- ' u -0.50 

-5.15 ±0.80 

-7 10 + ' 40 
' lu -0.35 

-5.30 ±0.20 


-4.0 

-3.16 

-4.46 

-5.6 

-4.45 


S 


14.80 ± 0.20 


-5.30 ±0.20 


-4.8 


Mn 


12.60 ± 0.10 


-7.50 ±0.15 


-6.58 


Fe 


14.10 ± 0.15 


-6.00 ±0.15 


-4.50 


Ni 


12.65 ±0.30 


-7.45 ±0.30 


-5.72 


Zn 


12.40 ± 0.30 


-7.70 ±0.30 


-7.4 



Table 3. Hydrogen and deuterium equivalent widths. Wave- 



lengths and oscillator strengths are taken from Morton ( 1991 ). 
H I and D I Ly e have been corrected for the influence of Fe II 
absorption at 937.652 A, Ff I and D I Ly ( for the influence of 
the H 2 We Q(l), 4-0 line at 930.574 A 



Ion 


Lyman line 


Wavelength [A] 


log /A 


W x [mA] 


Hi 


S 


949.7431 


1.122 


337.0 ± 20.1 




e 


937.8035 


0.864 


218.0 ± 21.0 




c 


930.7483 


0.651 


240.3 ± 18.9 




■'/ 


926.2257 


0.470 


222.8 ± 17.0 




K 


919.3514 


0.043 


202.1 ±21.5 


Dl 


s 


949.4847 


1.122 


84.6 ± 10.0 




e 


937.5484 


0.865 


47.2 ± 17.5 




c 


930.4951 


0.651 


24.8 ± 16.1 




■'/ 


925.9737 


0.470 


24.4 ± 11.7 




K 


919.101* 


0.043 


< 30.7 



Wavelength calculated from H I Ly k. 



Transition A [A] 


r 


Wx 


mA 


Notes 6 


J = 0( 5j = l,£j = 0);log 


Nj = 


15.10 ±0.20 





Ly R(0), 0-0 


1108.128 


0.00173 


18.6 ±4.8 




Ly R(0), 1-0 


1092.194 


0.00596 


27.7 ±4.0 




Ly R(0), 2-0 


1077.138 


0.0119 


46.7 ±6.1 




Ly R(0), 3-0 


1062.883 


0.0182 


41.1 ±6.0 




Ly R(0), 8-0 


1001.826 


0.0266 


< 34.0 


b 



J = 1 (9 J = 9,£j = 0.01469); log Nj = 15.50 



+0.20 

0.30 



We Q(l), 0-0 


1009.772 


0.0238 


53.2 ±7.1 




Ly R(l), 1-0 


1092.732 


0.00403 


22.6 ± 5.1 




WeQ(l), 1-0 


986.798 


0.0364 


44.7 ± 7.9 




Ly R(l), 2-0 


1077.698 


0.00809 


53.2 ± 7.2 




WeQ(l), 2-0 


966.097 


0.0349 


< 56.8 


b 


Ly P(l), 3-0 


1064.606 


0.00584 


15.6 ±5.0 




Ly R(l), 4-0 


1049.958 


0.0160 


36.1 ±4.9 




Ly P(l), 7-0 


1013.434 


0.0205 


< 40.2 


b 



J = 2 (gj = 5,Ej = 0.04394); logiVj = 14.10 



+0.15 
0.20 



We Q(2), 0-0 
Ly R(2), 1-0 
We Q(2), 2-0 
Ly R(2), 4-0 
Ly R(2), 5-0 


1010.941 
1094.244 
967.278 
1051.497 
1038.690 


0.0238 

0.00367 

0.0348 

0.0147 

0.0170 


22.3 ±6.1 
< 10.6 
15.5 ±5.5 n 
20.7 ±4.6 
9.7 ±3.8 


J = 3 (9 j = 


21, Ej = 0.08747); logA^ 


j = 13.95±2;» 


Ly P(3), 1-0 
Ly P(3), 3-0 
We Q(3), 4-0 
Ly P(3), 5-0 


1099.788 
1070.142 
933.581 
1043.498 


0.00240 
0.00721 
0.0193 
0.0104 


< 5.2 n 
13.2 ±4.3 
< 13.7 
11.8 ±3.5 


J = 4 (gj = 


9,Ej = 0.14491); logiVj 


= 13.80i-° : » 


Ly P(4), 3-0 
Ly Q(4), 4-0 
Ly P(4), 6-0 


1074.313 
935.959 
1035.184 


0.00733 

0.0192 

0.0108 


12.2 ±3.6 
9.4 ±5.6 
< 11.4 


J = 5 (gj = 


33, Ej = 0.21575); log AT/ < 13.40 


We Q(5), 1-0 
We Q(5), 2-0 
Ly P(5), 5-0 


994.924 
974.286 
1052.499 


0.0361 
0.0346 
0.0105 


< 9.0 
< 10.7 

< 3.9 


J = 6 (9 j = 


13, Ej = 0.29940); log A 


j < 13.40 


Ly R(6), 1-0 
Ly P(6), 4-0 
Ly R(6), 6-0 
Ly P(6), 6-0 
Ly R(6), 8-0 


1109.859 
1071.497 
1041.711 
1045.805 
1019.101 


0.00307 

0.0102 

0.0143 

0.0122 

0.0133 


< 4.1 

< 6.8 

< 4.8 

< 5.9 

< 3.7 


J = 7 (gj = 


45, Ej = 0.39521); log A 


j < 13.50 


Ly R(7), 5-0 
Ly P(7), 6-0 


1059.992 
1052.386 


0.0134 
0.0123 


< 7.7 

< 4.9 



a f values from Morton & Dinerstein ( 1976 ) 

b b: line is possibly blended; n: spectrum has locally small S/N 



